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Abstract. The question, what is the role of freeze-out of chemical species 
in determining the molecular abundances in the interstellar gas is a matter of 
debate. We investigate a theoretical case of a dense interstellar molecular cloud 
core by time-dependent modelling of chemical kinetics, where grain surface reac- 
tions deliberately are not included. That means, the gas-phase and solid-phase 
abundances are influenced by gas reactions, accretion on grains, and desorp- 
tion, only. We compare the results to a reference model, where no accretion 
occurs and only gas-phase reactions are included. We can trace that the purely 
physical processes of molecule accretion and desorption have major chemical 
consequences on the gas-phase chemistry. The main effect of introduction of 
the gas-grain interaction is long-term molecule abundance changes that come 
nowhere near an equilibrium in during the typical lifetime of a prestellar core. 
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1. INTRODUCTION 

The observed abundances of various species in the interstellar medium are 
frequently explained by gas-phase reactions. However, in dense and dark clouds 
elements other than H and He mostly accrete on dust-grains, forming icy mantles. 
There they may undergo chemical reactions and photoprocessing, while still in an 
adsorption-desorption equilibrium with the gas phase. The extent and significance 
of the surface reactions are matter of debate (Woodall et al. 2006). Models 
that include grain surface-phase chemistry in addition to gas-phase chemistry e.g. 
Hasegawa et al. (1992), Allen & Robinson (1977), generally have better and wider 
possibilities to reproduce the production of various species. However, the solid 
phase itself acts as a reservoir, which withdraws molecules from active gas phase 
chemical transformations. We aim to investigate the effect of the physical freeze- 
out of molecules onto grains on gas-phase chemical processes. A theoretical case is 
examined by chemical kinetics modelling, where the freeze-out of heavy elements 
occurs but the surface reactions are not included, with H2 production on grains as 
the only exception. Therefore, in addition to the gas-phase reactions the accretion 
and desorption (Sect. 2.1) of species onto and from the grains are taken into 
account, only. 
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Table 1. The adopted abundances and depletion of elements used in the model. 





Total 






Initial 


Element 


abundance 3. 


Depletion, 


Source 


gas-phase 






% 




abundance 3 


H 


1.0 







1.0 


He 


0.1 







0.1 


C 


2.9E-04 


39 


J09 c 


1.8E-04 


N 


7.9E-05 


22 


J09 


6.2E-05 


O 


5.8E-04 


42 


J09 


3.3E-04 


Mg 


4.2E-05 


95 


J09 


2.3E-06 


Fe 


3.5E-05 


99 


J09 


2.0E-07 


S 


2.1E-05 


92 


J09, P97 d 


1.8E-06 


Na 


2.1E-06 


75 


P97, PPG84° 


5.3E-07 



a relative to hydrogen 

b Prodanovic et al. (2010) 

c Jenkins (2009), with F* = 1 

d Pagel (1997) 

c Phillips et al. (1984) 

2. THE MODEL 

To investigate the effect of an inert reservior on molecular abundances we compare 
the results of two chemical kinetic models: Model 1 with gas phase species taken 
into account only, and Model 2, where some of the neutral species are assumed 
refractory and can stick to grain surfaces, and the frozen grain mantle forms. The 
model is essentially a time-dependent and simplified version of the model described 
by Kalvans & Shmeld (2010) and Kalvans & Shmeld (2012) - a plane-parralel in- 
terstellar cloud, illuminated from both sides by the interstellar ultraviolet and 
cosmic-ray radiation field. However, we examine only a deep cloud core, where 
the visual absorption is 20 magnitudes, i.e. the interstellar ultraviolet radiation 
plays essentially no role. We use the UDFA06 dipole (Woodall et al. 2006) ki- 
netic chemical reaction set. The gas temperature is taken 15K, dust temperature 
10K. We assume atomic initial abundances. The calculations are performed with 
S.OtimeslO 11 s as the initial time in 15 time steps with 1.4 as the multiplying factor 
until 3.3iimesl0 13 s is reached as the final time. For Model 2 accretion on dust 
grains and desorption of molecules but no surface or other solid-phase reactions are 
included. Of the 301 chemical species included 83 are assumed refractory and form 
a mantle. The abundance and depletion of the elements for the model is shown in 
Table 1. We opted for a case, when elements are already heavily depleted, since 
the case of an already formed molecular cloud core is investigated. 

2.1. Accretion and desorption 

The rate coefficients for accretion of molecules onto grains (Willacy & Williams 
1993), thermal evaporation and desorption by grain heating by cosmic rays (Hasegawa 
& Herbst 1993), desorption by cosmic-ray induced photons (Prasad & Tarafdar 



The effect of inert molecule reservoir on interstellar abundances 



279 



1983, Willacy, Williams 1993), desorption of surface molecules by energy released 
by H2 molecule formation on grains (Roberts et al. 2007), and direct ejection of 
molecules from the mantle by cosmic-rays all are calculated according to Kalvans 
& Shmeld (2010). The species desorbed by H2 formation are taken those with 
binding energies E\, < 1210-ftT (binding energies from Aikawa et al. 1997). 

2.2. The formation of the mantle 

The mantle-surface model is rather simple and was made to include the two phases 
in a way that in the period modelled they essentially accrete with only minor return 
to the gas phase. The mathematical model is such that the surface molecules 
undergo various desorption processes (Sect. 2.1) and turn into mantle molecules 
with a constant rate, while the mantle molecules may undergo direct ejection into 
gas phase and return to surface. First, the rate coefficient of mantle molecules is 
calculated by 

^dig—up I'cr X f , (1) 

where t cr = 10 12 (Bringa & Johnson 2004) is the average time between two 
successive strikes by Fe cosmic-ray nuclei and v is the assumed number of strikes 
needed to return all mantle molecules to the outer surface for some time. Then, the 
rate coefficient of surface molecules turning into mantle molecules is assumed to be 
kmantie = 100kdig-up and it would produce a mantle-to-surface total abundance 
ratio of 100:1. 

We choose w=5000, so that the time to reach the approximate 100:1 ratio is on 
the order of 1 Myr, more or less consistent with the time of mantle compression, 
as found by Palumbo (2006) and Accolla et al. (2011). A high abundance of 
molecules exposed to outer surface (low mantle-to-surface ratio) means a high 
porosity of the mantle. The porosity is reduced over time by cosmic-ray hits, 
exothermic reactions and other events (see Palumbo, 2006, Accolla et al., 2011, 
Raut et al., 2007a, b, 2008). 

This mathematically simple picture is the same as used in Kalvans & Shmeld 
(2012). It is not quite physically adequate, and is a little more sophisticated than 
the one used in Kalvans & Shmeld (2010). However, it includes the accretion pro- 
cess and interaction of the solid grain mantle with gas phase and allows us to study 
the effect of the introduction of an inert reservior on the gas-phase abundances of 
chemical species. 

3.. RESULTS AND DISCUSSION 

The results are graphically presented in Figs. 1 to 6. The calculated abun- 
dances for selected species for Model 1 are shown in Fig. 1, the gas, solid and 
total abundances for Model 2 are shown in Figs. 2, 3 and 4, respectively. For 
comparison, the ratio for Model 2 to Model 1 calculated gas and total abundances 
are shown in Figs. 5 and 6, respectively. 

The model with constant physical parameters approximately represents a case 
of a slow-contracting molecular core. Chemical reactions compete with the ac- 
cretion of molecules on grains. Small molecules form and freeze faster than more 
complicated and heavy species. However, they have smaller adsorbtion energies 
and more easily desorb to the gas-phase. At time 1 Myr > 95 % of the heavy 
elements are in solid phase and the freeze-out continues, consistent with models 
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Fig. 1. Calculated gas-phase relative abundances for selected species for Model 1 (pure 
gas- phase). Pay attention to the relative abundance (ordinate) scale for each graph. 



by Brown & Charnley (1990), Hasegawa et al. (1992) and Roberts et al. (2007). 

The molecule accretion process facilitates the transformation of free atoms into 
molecules. In Model 1 atoms like C, O, N and S are the major species along with 
the CO molecule (Fig. 1), while in Model 2 molecules like CO, 2 , H 2 0, CN, 
NO, S 2 etc. take over, especially so, when one takes into account the total (gas + 
solid phase) abundances (Figs. 2, 3 and 4). Simple molecules with usually 2 metal 
atoms contain the majority of heavy elements. The largest increase of abundances 
for Model 2 compared to Model 1 at 1 Myr is for molecules with 1-3 metal atoms 
with varying degrees of hydrogen content. The abundance increase for multi-atom 
carbon-chain species is usually only moderate and they have become depleted 
from the gas-phase. These results make it difficult to explain high abundance of 
these molecules in quiescent clouds without grain reactions (as in Herbst & Leung 
1989). As the accretion-desorption equilibrium approaches, the total abundance of 
many important small molecules (CO, CN, N 2 , CS, etc.) begins to slowly become 
smaller. The atoms are then distributed to species with longer formation and 
freeze timescales, including multi-atom species. Thus, the chemical diversity for 
long lived clouds can be enhanced with time by this process. 

Except for a few molecules (N 2 , NH 2 CN) the change of concentration for Model 
1 molecules has nearly ceased already at 0.5 Myr, indicating a state close to a 
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Fig. 2. Calculated gas- phase relative abundances for selected species for Model 2 
(with freeze processes). Pay attention to the relative abundance (ordinate) scale for each 
graph. 



chemical equilibrium. It is very unlike Model 2, where both, gas and solid phase 
abundances are still considerably changing at 1 Myr. One can conclude that the 
physical gas-grain interaction can cause lasting chemical effects that may allow 
to deduce the history of an observed molecular cloud core or young stellar ob- 
ject. However, for a credible investigation a full interstellar gas-grain physics and 
chemistry model is needed. 
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Fig. 5. The comparison of calculated gas- phase abundances for selected species be- 
tween Model 2 (with freeze processes) and Model 1. The ratio of Model 2 to Model 1 
gas abundances is shown. Pay attention to the M2/M1 (ordinate) scale for each graph. 
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Fig. 6. The comparison of calculated total (gas + solid phase) abundances for selected 
species between Model 2 (with freeze processes) and Model 1. The ratio of Model 2 to 
Model 1 total abundances is shown. Pay attention to the M2/M1 (ordinate) scale for 
each graph. 
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